Penetration of hot electrons through a cold disordered wire. 
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We study a penetration of an electron with high energy E 3> T through strongly disordered wire 
of length L 2> a (a being the localization length). Such an electron can loose, but not gain the 
energy, when hopping from one localized state to another. We have found a distribution function 
for the transmission coefficient T. The typical T remains exponentially small in L/a, but with the 
decrement, reduced compared to the case of direct elastic tunnelling: InT « 0.237 ■ 2L/a. The 
distribution function has a relatively strong tail in the domain of anomalously high T; the average 
T oc (a/L) 2 is controlled by rare configurations of disorder, corresponding to this tail. 

PACS numbers: 72.20.Ee, 73.21.Hb, 73.63.Nm 



The electronic transport in disordered one-dimensional 
systems was extensively studied in the past 50 years [2, El]. 
In the non-interacting system all the states are localized 
Q , so that the transmission coefficient T of a finite sys- 
tem exponentially decays with L. T varies from sample 
to sample, since it depends on the configuration of dis- 
order. For strongly disordered chains the distribution of 
s = —(l/a) InT is a narrow gaussian: 

^o(s)ocexp{-( S -l) 2 /(A,s) 2 }, ^s = Ba- 1 ' 2 , (1) 
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FIG. 1: The setup. A disordered wire of length L is in a good 
contact with two reservoirs, the entire system is in equilib- 
rium. A small flux I[ n of electrons with energy E S> T is sent 
to the left end of the wire through the left reservoir, and the 
part Ir = Tliu finally reaches the right reservoir. 



a = 2L/a > 1, 



(2) 



being the principal large parameter of the theory (the 
measure of the localization strength), B ~ 1 being the 
model-dependent factor [J]. Thus, the conductivity of a 
noninteracting system is zero at L — >• oc. 

The interactions (e.g., with phonons) lead to a finite 
equilibrium conductivity a of the hopping type (see [HJ). 
At some L ~ L (T) the exponential L-dependence of the 
conductance G is changed to G — o~/L, with a, expo- 
nentially dependent on the temperature T of the system. 
The specifics of strongly disordered Id systems was prop- 
erly taken into account in [g, 0] ; it was shown that the 
conductivity in the variable range hopping regime is con- 
trolled by rare fluctuations of the density of states at the 
fermi level - the "breaks" ; as a result 



cr « exp{-T /2T}, T = 1/ga, 



(3) 



where g is the average density of states. The result (|3J) 



c (T /T)VC«H-i) }j 
), where electrons 



does not obey the Mott law a oc expj 
valid in dimensions d > 2 (see jj 
can easily circumvent the breaks. The transport in the 
strong field £ was studied in 

I, EMI- Here the current 

I oc exp{— 8To/e£ a} strongly depends on £ , not on T. It 
is insensitive to the breaks, but, on the other hand, the 
distribution of electrons is far from equilibrium. 

In the present paper we study a different situation, 
where the current through the system arises due to a 



small group of strongly nonequilibrium high-energy par- 
ticles, so that the occupation numbers of most electronic 
states remain essentially in equilibrium. 

A disordered wire of length L is in equilibrium with two 
metallic leads (reservoirs) at temperature T and chemical 
potential ep (see FigJTJ). In the left reservoir, however, 
a small amount of nonequilibrium particles with energies 
E ^> T (E is measured with respect to ep) is injected, so 
that the current I lVL reaches the left end of the wire. 

Hot electrons, injected into the wire, weakly interact 
with the thermal bath, their energy is not conserved. 
However, as long as E ^> T, only the processes in which 
the energy is transferred from electron to the bath, not 
vice versa, are allowed. This is only true for not very 
long wires L < L Mott {T), where L Mott (T) ~ a(T /T) 1 / 2 
is the length of the typical Mott hop at given T. Under 
this condition the equilibration does not have chance to 
occur before the electrons escape from the wire. In this 
letter we also do not take into account correlation effects 
(like Coulomb gap) due to electron-electron interactions. 

Obviously, only the localized states with the energies 
Ej in the interval between the Fermi energy and the ini- 
tial energy E of the injected electron are relevant for our 
problem. We enumerate them according to their energies: 
< E\ < £2 <...<£„< E. These "quasiresonant" lev- 
els play an important role in the transport physics, as 
the electron, travelling through the chain with initial en- 
ergy E, can make intermediate stops only at these sites. 
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Indeed, all the sites with e; < are occupied, while the 
sites with a > E cannot be reached, as no energy can 
be absorbed. The spatial positions of quasiresonances 
are ^ = Lxi, the independent random variables Xi are 
homogeneously distributed in an interval < Xi < 1. 
The number n of quasiresonances is itself a random vari- 
able, described by the poissonian distribution p{n 1 N) — 
N n exp(— N)/nl , where N = LEg is the average number 
of quasiresonances. Thus, each wire is characterized by 
"the configuration" C = {n, {x\, ...!„}} and the average 
of any C-dependent quantity A(C) over the ensemble of 
wires is A = ^2^L Q p(n, N)Y\^ =1 dxtA(C). In partic- 
ular, the distribution function P/v(s) = S(s — s(C)). In 
this letter we focus on the most interesting case N 3> 1, 
when the poissonian distribution is sharp and one can 
simply average over at fixed n w N . 

How can an electron get from left reservoir to the right 
one? Besides the obvious possibility of the direct elastic 
tunnelling (Figj2^,), there are also numerous "inelastic 
staircases" (FigJ2b,c), in which an electrons makes in- 
termediate stops at certain quasiresonant states, while 
the excess energy at each hop is transferred to the ther- 
mostat. Each staircase S is characterized by the choice 
of a subset of K (0 < K < n) quasiresonances, with 
e kl < s k2 < ■ ■ ■ < e kK and x kl > x k . 2 > ■ ■ ■ > x kli . 

Each staircase contributes to the transmission: 



T{C) = T{S\C) oc e 
s 



s(C) 



T(S\C) cx e 



-as(S\C) 



(4) 



where the summation runs over all the staircases, pos- 
sible for given configuration C. Under the condition @ 
the sum in @ is dominated by only one - the optimal 
- staircase <S op t(C) that corresponds to minimal s(S\C), 
so that s(C) ~ s[S op t{C)] — mmg s(S\C). In a typical 
situation the longest hop in the optimal staircase is the 
last one, then goes the last but one, etc. Therefore the 
value of s(S\C) is controlled by few last hops in S, while 
the multitude of short hops in the upper part of the stair- 
case are of only secondary importance. The most natural 
assumption about the structure of the optimal staircase 
would be the scaling hypothesis: the distribution func- 
tion Pfe(^fc) for random variables £ k — x ik /x ik l does not 
depend on k. Such a simple self-similar structure was, 
however, not observed in our numerical experiments: V k 
manifestly depended on k. 

The explicit expression for s(S\C) can be found with 
the help of the stationary master equation for the 
quasiresonant levels populations fa: 



0, 



(5) 



The first term in ([5|) is the incoming flux of particles from 
the left reservoir; the second term describes the particles, 
coming to the level i from all other levels with higher 
energies (hence the condition j > i); finally, the third 
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FIG. 2: (color online) Possible scenarios for electron passing 
from the left to the right reservoir. "Natural hops" (occurring 
with the probability close to unity) are shown by thick green 
arrows, the non-natural hops (that happen with exponentially 
small probability) are shown by thin red arrows, (a): Direct 
elastic tunnelling; (b,c): Inelastic staircases with intermediate 
stops at quasiresonances. 



term takes into account all possible escapes from the z-th 
level: P^out = Pi^L + Pi^R- 
transmittance is 



V ,. , /', . ,. The system's 



T = I R /I- m = Pl 



(6) 



The rate of transitions between j — > i is Pj^i — 
Pj^iisj — Ei)e~ a \ Xi ~ Xj ' . Matrix elements of the electron- 
thermostat interaction, entering fi^j are smooth power- 
law functions of the energy transfer Ej — £j. It means, 
that if we are interested only in the exponential depen- 
dencies, we do not have to take these matrix elements 
into account. Thus, in the exponential approximation 
we can write 



Pi^ j( xe{i-j)e- a \'*-^, Pi 



— >out c< e 



(7) 



Xi{C) = mini x tl 1 - a*, min{|a;j - xj\} > , (8) 



being the distance from the i-th quasiresonance to its 
"natural descendant" - a closest neighbor with lower en- 
ergy, or to one of the two reservoirs. The solution of 
the system of equations ^ can be written in a recurrent 
form: 

fi — {Pl— ril Pi— >out)^in + ^ ' fj{Pj^i/ Pi— ►out ), (9) 

j>i 

which allows for finding /, provided all fj with j > i 
are already found. According to the exponential approx- 
imation, justified by the large parameter ([2]), any sum, 
occurring in ([6]) or in (|9]), is dominated by a single term 
with the smallest negative exponent. As a result, the 
normalized probability for the electron to make a hop 
i -t j is p(i -t j) = P^/P^out = exp{-aAs(i -> j)}, 
where As(i — > j) = \xi — xj\ — Xi- Then, having in mind 
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FIG. 3: (color online) A "fortunate" configuration (with m = 
5). The last m hops in an optimal staircase are "natural". 
There are no quasiresonances in the shaded domain. 



that T(S\C) — T[k=i P(ik — > ifc-i)) we arrive at 

K+l 

s(S\C) = ^ As(i k -> ik-i), 



(10) 



fc=i 



so that As(ik — > ik—i) has the meaning of the contribu- 
tion of the fc-th hop to the exponent s. In (TTOl) xq = 1 
and Xif+i = are the positions of the right and the 
left reservoirs, correspondingly; xk+i — 0. From (fTU|) 
it immediately follows that s(S\C) = 1 — 



K+l 
k=l 



Xi 
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principle, our problem can be solved by means of enumer- 
ation of all staircases, possible for a given configuration, 
and choosing the optimal one. The direct solution of the 
master equations, however, leads to the same result. 

Note, that As(ik — > ik-i) = if the ik-i-th quasires- 
onance is the natural descendant of the i^-th one. The 
corresponding hops we will call "natural hops" in what 
follows. Clearly, to minimize s, it would be nice to have 
a staircase, where all the hops iu — > ik-i ( or & t least as 
many of them, as possible) are natural. We will see, how- 
ever, that such a "natural staircase" is possible to find 
only for some rare "fortunate configurations" . 

Suppose that certain configuration C generates an op- 
timal staircase, in which a sufficiently long subsequence 
of last to hops k m —> k m —i fci — > R is "nat- 

ural", see Figj3] It means that, if an electron somehow 
manages to get to the upper level k m in this subsequence, 
then it makes its descending way to the right reservoir 
through the rest of the staircase with a probability that 
is close to unity. Then, if this upper level is close enough 
to the left reservoir (namely, if Xk m < s), then one can 
guarantee that s(C) < s, no mater whether the preced- 
ing short hops in the staircase are natural or not. The 
necessary conditions for a given configuration C to be a 
fortunate one are as follows: The last jump (to R) should 
start from the lowest level (i.e., k\ = 1), and this level 
should be localized in the right half of the wire: x\ > 1/2. 
In general, the level k p in the optimal staircase should 
be the lowest one among all the quasiresonances with 
< Xi < Xk p _ 1 ; this level should be in the right half of 
the stretch: Xk p _ 1 /2 < Xk p < Xk p _ 1 - This should be true 
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FIG. 4: (color online) (a): The "natural path" consisting of 
hops to nearest neighbors with lower energy. The particle 
starts from the middle of the wire, so that in one half of the 
ensemble of samples the natural path leads the left reservoir, 
and it leads to the right one in the other half, (b): The 
particle is initially placed near the left reservoir, the natural 
path leads there in almost all samples. To get to the right 
reservoir, a particle has to make non-natural hops. 



for all jumps with p < m, where m is determined by the 
condition Xk m < s. 

The probability for all the m last jumps in the opti- 
mal staircase to be natural is 2~ m . However, m is not 
determined solely by s: for fixed s m still fluctuates from 
configuration to configuration. 

Let us introduce statistically independent random vari- 
ables £ p = 2[x kp - x kp _ 1 ]/x kp _ 1 , homogeneously dis- 
tributed in the interval < £ p < 1. Then s = n^Li(l + 
£ P )/2. The distribution P m {<t>) of the random variable 
<t> = to In 2 — ln(l/s) = m (l + £p) can ^ e obtained 

with the help of the Fourier transformation: 



P m (4>) 



[ * e -** 


- 2*9+1 _ 1 - 


rn 


/ 2tt 


iq + 1 





Fcois) w / d<t> p m{<f>)S(s - exp{0 - to In 2}) = 

_ P m (mln2-ln(l/s)) _ f dq (1 + iq)e iqla{ - 1 /^ 



s2 ri 



2vrs 2-( 1 +'"?) + iq 



where b ~ 1 is a universal constant. Although this 
asymptotics describes only a small fraction of the con- 
figurations, it turns out to be sufficient for finding the 
average transmission coefficient: 



Tt 



bs exp{— as}ds — ba 2 oc L 



(11) 



The hopping motion of hot particles, accompanied by 
the emission of energy, earlier was studied in the context 
of recombination of photo-excited electron- hole pairs [l3[ . 
For that end it was sufficient to take into account only 
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N=1200 




FIG. 5: Distribution functions of the logarithm of trans- 
parency for different values of average number of quasires- 
onaces N = ELg. For N > 200 the curves practically 
merge and coincide with Foo(s), shown by thick line. The 
^-dependence of the average logarithm of transmittance is 
shown on the inset, the solid line being the asymptotic ex- 
pression (fT2")l . 



the natural hops since a particle was assumed to be cre- 
ated in the bulk, far from the ends of the sample, so that 
it could not escape to a nearby lead (see FigHk). Under 
these conditions the evolution of particle density profile, 
averaged over the ensemble of samples could be described 
by a peculiar diffusion, in which the distribution of hops 
lengths is rescaled by a fixed parameter q > 1 after each 
hop. If one would extend the same approach to the setup 
where the particle is initially placed near one of the ab- 
sorbing leads, then one easily finds for the probability to 
reach the opposite lead (and thus the average transmis- 
sion coefficient) T oc L - * 3 , where the exponent (3 (as well 
as the parameter q) is universal. Note that this result 
is in agreement with (|11[> . However, as we have already 
mentioned, T is controlled by very rare anomalous sam- 
ples with high transparency, while in a typical sample 
the natural path would lead to the nearby left lead (see 
Fig lib)- Therefore, to find the probability to reach the 
right lead in a typical sample one should take into ac- 
count the non-natural hops, that were ignored in flU . 
The "diffusional approach" 



13j, being an adequate in- 



strument for finding 7", is useless for the determination 
of the distribution of T. 

The distribution Fn(s) for general s and N can only be 
found numerically. We generated an ensemble of ~ 10 6 
random configurations C and calculated corresponding 
s(C) with the help of the recurrent formula ([9]). The 
results of our Monte-Carlo simulations are summarized 
in Fig[5] The distribution functions F N (s) are wide: the 
dispersion of s is of order of s for all N. For N > 300 
practically Fn(s) s» F oq (s), and ^(s) indeed shows a 
linear low-s asymptotics with b ks 29. 

The average value s(N) monotonically decreases with 
N, tending to a finite limit s(co) s» 0.237 at N — >■ oc. The 
convergence is, however, extremely slow: for N = 1000 
the relative difference is still of order of 2%. The conver- 



gence can be improved dramatically if one introduces an 
asymptotic correction according to the empiric formula 

s(N) « 0.237 + 0.598 In N/N, for N > I . (12) 

The deviation of experimental s(N) from the asymp- 
totics ([nj is less than 2% already for N ~ 30. The N- 
dependent corrections are controlled by many short hops 
in the beginning of the optimal staircase. The scaling hy- 
pothesis, mentioned above, would lead to s(N) — s(oo) oc 
1/N, which is inconsistent with our numerical data. This 
is another argument against the scaling. A fit of our nu- 
merical data for the -ZV-dependence of an average number 
K of hops in the optimal staircase gives 



K{N) w 0.39(ln7V) 2 + 2.4, 



for N » I 



(13) 



which is perfectly consistent with the result (|T2l and, 
again, inconsistent with the scaling hypothesis, (the lat- 
ter would give K(N) oc IniV). The elucidation of the 
structure of the initial part of the optimal staircase and 
the origin of the empirical laws (|f 211 31) still remains a 
challenge. 

In conclusion, we have found the distribution function 
of the transmission coefficient for the inelastic penetra- 
tion of a cold disordered wire by a hot electron. The 
applications of these results to specific physical effects 
will be presented in a long paper to follow. 

We are indebted to M.E.Raikh for valuable comments. 
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